Surface superconducting states in a polycrystalline MgB2 sample 
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We report results of dc magnetic and ac linear low-frequency study of a polycrystalline MgB2 
sample. AC susceptibility measurements at low frequencies, performed under dc fields parallel to 
the sample surface, provide a clear evidence for surface superconducting states in MgB2. 

PACS numbers; 74.25.Nf, 74.25.0p, 74.70.Ad 



Nucleation of superconducting (SC) phase in the sur- 
face sheath under a dc magnetic field parallel to the sur- 
face was predicted by Saint-James and de Gennes more 
than 40 years ago 1\. This prediction was made in the 
frame of the one-band isotropic Ginzburg-Landau (GL) 
model and experimental confirmations of this prediction 
were described in various publications (see, for exam- 
ple, The discovery of superconductivity in two-band 
anisotropic MgB2 raised a question about the existence 
of surface superconducting states (SSS) in this material. 
Transport measurements under dc fields indicate that the 
onset of the SC transition occurs above the second criti- 
cal field, Hc2 determined from specific heat data QQ. 

In this work we present the experimental results of dc 
and ac measurements on polycrystalline MgB2 sample. 
It appears that large losses with a maximum at certain 
temperature dependent dc field and complete diamag- 
netic screening of an ac field are observed in magnetic 
fields when the dc magnetic moment is very small. This 
is clear evidence for the existence of SSS in this material. 

The MgBa sample with ^ 37.5 K and AT^ = 0.5 K 
was prepared using the Hot Isostatic Pressing method . 
The sample dimensions are 10 x 3 x 1 mm'^. DC magne- 
tization curves were measured by a SQUID magnetome- 
ter, ac susceptibility was measured by the pick-up coils 
method in the frequency range 5 < lo/2ti < 1065 Hz. 
An ac "home-made" setup was adapted to a commercial 
SQUID magnetometer and its block diagram was pub- 
lished elsewhere 0- The experiments were carried out 
as follows. The sample was cooled down in zero mag- 
netic field, then both dc (Hq) and ac magnetic field with 
amplitude ho were applied and the amplitude and phase 
of ac response was measured. Both Hq and were par- 
allel to the longest sample axis. 

Assuming that in zero dc field there is a complete 
screening of ac field without any losses, one could mea- 
sure the absolute value of ac susceptibility, x' and x"j 
as a function of external parameters, such as dc field, 
frequency, and temperature. At low ac amplitudes the 
response is linear and does not depend on ho as shown in 
the inset (a) to Fig.Ht, where the x" at w/27r = 1000 Hz, 
Hq = 4.3 kOe, and T = 36 K is plotted as a func- 
tion of ho. Fig. ^ shows magnetization loop M{Hq) 
at T = 34 K (upper inset) and the high field M(ffo) 



data in an extended scale (main panel), where the irre- 
versibility field, Hirr and Hc2 are indicated. The lower 
inset of Fig. ^ presents Hirr{T) curve. These results 
are consistent with the data reported previously (see, for 
example, (9| and references therein). 

The inset (b) in Fig. ^ shows dc and ac susceptibility 
measured at 34 K as a function of dc field. It is evident 
that the Xdc = M{Hq)/Hq is already very small at 5 kOe, 
whereas the x' becomes zero only above 12 kOe. The dc 
field dependence of x' and of x" is shown in Fig.^ (main 
panel). Here again, the x' curve exhibits ideal diamag- 
netism in fields much higher than 5 kOe. Note the high 
peak of losses at about 11 kOe, Fig.^, which consider- 
ably exceeds the losses both in the normal and Meissner 
states. A frequency dispersion of x is also evident (see 
Fig. |2Jl . All these features are typical for the ac response 
in SSS of low-temperature superconductors (LTS) 2, ^. 
However, there is an apparent difference. In isotropic 
LTS these features are observed only at Hq > Hc2 , where 
bulk dc magnetization is zero. In contrast, in polycrys- 
talline MgB2 SSS are observed below Hc2 (see Figs. ^ 
and ^p) and coexist with dc bulk magnetization. The 
dc moment in this region is reversible and small, but un- 
ambiguously has a bulk origin. In order to prove this, 
we applied a low-frequency shaking field with the am- 
plitude of 2.5 Oe. We did not observe any change of 
the dc magnetization curve. Any possible nonzero sur- 
face current contribution to dc magnetic moment has a 
nonequilibrium character if the volume is in the normal 
state, and one could expect hysteresis and disappearance 
of the dc magnetization after the application of the shak- 
ing ac field. On the other hand, the ac response cannot be 
attributed to bulk properties. For example, at T = 34 K 
and Hq = 10.2 kOe x' at 1065 Hz is about 0.5 (the mid- 
point of the transition in Fig.^, main panel). This value 
is 3 orders of magnitude higher than the experimental 
dM/dHo measured under the same conditions. 

The frequency dependence of x' and x" at constant dc 
fields is shown in Fig. |21 It is clear that in some fields 
x'{^) oc ln(ti;). x" show a maximum at Lojl-K = Fm and 
this maximum moves with dc field as shown in Fig. |31 
One may predict that, similar to the spin-glass state, 
the low frequency dynamics of SSS could be character- 
ized by a broad spectrum of relaxation times with some 
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FIG. 1: (Color online) (a) Magnetization curve, M{Ho), in 
high magnetic fields. V - theoretical calculations for two-band 
model. Upper inset - magnetization loop after ZFC at T = 
34 K. Lower inset - HirriT). (b) Magnetic field dependence of 
x" and x' ^-t <^/27r = 1065 Hz in high magnetic fields. Insets: 
(a) ^'{ho) dependence at T = 36 K; (b) Xdc = M{Ho)/Ho 
for the virgin curve, and x'iHo) in whole range of dc field. 
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FIG. 2: (Color online) Frequency dependence of x" (a-) a-nd 
x' (b) at T = 34 K for different magnetic fields. The symbols 
for the fields in both figures are identical. 
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maximum value r. Indeed, Fig.|31shows that Fm{H^) di- 
vides the frequency - magnetic field phase diagram into 
two parts, above and below the line, where ujt > 1 and 
wr < 1, respectively. 

GL equations for two-band anisotropic superconduc- 
tors were discussed in numbers of articles [lol ITlL IT^ 
and in linear approximation these take the form: 



FIG. 3: (Color online) Magnetic field dependence of Frr. 



(1) 



where the index s — 1(2) corresponds to the a{TT) bands, 
Hj = iVj + 2eAj/'hc, A - vector potential, As - order 
parameter of s band, p,s - inverse mass tensor and Sss' is 
different from zero only if s s' . The coefficients in these 
equations, Jls, cts{T) and 5ss' , are given through the BCS 
superconducting coupling matrix Xnm and could be found 
in 0, 01 . Let us consider the superconducting slab of 
thickness 21) in a magnetic field parallel to its surface and 



chose the coordinate system with a;-axis perpendicular 
to the surface, magnetic field along z-axis and the plane 
X = in the center of the slab. The second, Hc2, and the 
third critical magnetic fields, H^s, are determined as the 
maximum fields for which the eigen solutions of linearized 
Eq. ^ with vector potential A (0, Hx, 0) and appropri- 
ate boundary conditions could be found. Looking for the 
solution of Eq. |^ in the form: 



fs{x) exp{iky + ikzZ-ifisxzkzX + i^sxyhx^ /2) (2) 
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one obtains 



d^fs 



i^J-syy -^l^y)h'^{x - Xsffs + (« + g) f s + 

Sss' fs'expiicps') = 
g = „ _ni < 
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with 01 = -02 = (pixy - f^2xy){kx - hx'^/2) + {^i^^ - 

fJ-2xz)f^zX, and boundary conditions dfs{D)/dx — 0. Here 
the dimensionless variables are used: h = 2eHo/hc, 

f^sik — f^sik / f^sxx 7 — O^s / f^sxx : ^ss' — ^ss' / f-^sxx • ^c2 

determined by the maximum field for which the solution 
of Eqs. (O exists at fc = fc^ = with boundary conditions 
dfs{0)/dx = and dfs{D)/dx = and iJcs is determined 
by the solutions of this equation with k w hD, fs{0) = 
and dfs{D)/dx = 0. 

At first we discuss the one-band anisotropic super- 
conductor, 5ss' — 0. Equations ||2J) are of the same 
type as in the isotropic case with GL parameter k = 

^yy — Ji'^y, and using the known solution we obtain 

Hc3 — 1.695a(/i^a;/iyj^ — /i^j,)^^/^. Hc2 is also determined 

by Eq. ||2Jl and equals to a{fixxfJ'yy ~ fJ-xy 

)-i/2 J3]. Thus 

the ratio Hcj,/Hc2 for anisotropic crystals does not de- 
pend on the orientations of principal axes and is the same 
as for the isotropic case. For a two-band superconduc- 
tor with different /i^fc the analytical solution hardly could 
be found and the numerical simulations can be success- 
fully used. For an uniaxial crystal, i?c2 depends only on 
the angle 9 between the c-axis and the direction of the dc 
magnetic field. In polycrystalline samples only the grains 
with Hc2 > Hq give a contribution into the equilibrium 
magnetic moment .14t .15j : 



e{Ha2iO) - Ho)ei9)sin{e)de, 



(4) 



where function e{6) also takes into account the anisotropy 
of the vortex lattice and is equal to 1 for isotropic sam- 
ple, Pc = 1.16, and 0(j/) = 1 for y > and zero other- 
wise. For an uniaxial two-band superconductor, e(9) is 
not known yet while for the one-band superconductor it 
was calculated in |ll|, e{0) = [sin2(6') + 7^ cos^ie)]/-/^/^, 
where the anisotropy parameter 7 = j?r,2 (7r/2) /Hr,2 (0) ■ 
Using for the coupling matrix from [13 An = 0.81, 
A22 = 0.285, A12 = 0^19, A21 = 0.091 and /i2c/M2a = 
1.487, fiia/fi2a = 0.2 12] and varying fiic/nia, the equi- 
librium magnetization near Hc2 have been calculated for 
different temperatures. We found that using fiik for clean 
limit with fiic/lJ-ia — 0.048 and e — 1 gives the best fit 
to experimental data. The result of these calculations 
is shown in Fig. ^ by triangles. The GL parameters 
obtained from this procedure are: k = 18.6, 16.1 and 
13.8 for T = 34 K, 36 K and 37 K, respectively. The 
GL formulae for one band isotropic superconductor with 
K = Hc2/ Hc\/2, where He is the thermodynamic critical 
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FIG. 4: (Color online) The ratio Hc3/Hc2 as a function of 
polar angle 6 for <^ = 0, 30, 45, and 90 degree. Insets: (a) 
angular dependence of Hc2{0); (b) Hci,/ Hc2{f>) for 9 = 90°. 



field and M^q = {Hq - Hc2)/87r/?cK^, gives k = 7.7, 7.8, 
and 13.9 for the same temperatures. At T = 36 K and 
Hq — 3 kOe, Meg from this formula is 1.81 emu/cm'^, 
whereas the experimental value is significantly smaller, 
0.11 emu/cm^. The calculated Hc3/Hc2 ratio for a single 
crystal equals, with accuracy « 10%, a value of 1.69 pre- 
dicted in 1] for all orientations of principal axes, whereas 
both Hc2 and Hcs strongly depend on orientation. Fig.^ 
shows Hc3/Hc2 as a function of the polar, 0, and az- 
imuthal, ip, angles using the same set of material pa- 
rameters as above for T = 34 K that corresponds to 
T/Tc w 0.9. The inset (a) in Fig. |3| shows Hc2 as a func- 
tion of 9. The ratio Hcs/ Hc2{'fi) is shown in the inset (b) 
in Fig. ^ The calculated anisotropy parameter 7 = 2.9 
is in quantitative agreement with 7 reported in [T3I and 
somewhat larger than found in • 

The sample investigated in the present work is com- 
posed of randomly oriented grains with average size of 
about 40 nm and evenly distributed Hcs- The observed 
complete screening of ac field shows that the whole sam- 
ple surface is in SSS. The onset of the screening could 
be considered as the percolation transition and the ap- 
pearance of the large continuous clusters in SSS. Perco- 
lation transition in SSS of Nb was discussed in Ref . ■ 
Fig. [5I shows the expanded view of x' and x" at large 
fields and T = 34 K at two frequencies Ci;/27r = 105 and 
1065 Hz. The physical nature of the displacement of the 
curves along the _ff-axis with frequency is not clear yet, 
but the response at 1065 Hz gives a more accurate value 
of the field H* , at which the percolation transition takes 
place. The H* found from x' equals 11.5 kOc. The Hc2 
at T = 34 K is about 14.5 kOe. Using the calculated 
above 7=2.9, one could find that parallel to the c-axis 
is about 8.5 kOe and H*/H^^ « 1.35. Assuming 
random uniform distribution of c-axis orientation in the 
grains and Hc3{9ip) / Hc2{9) = 1.7 we calculated the frac- 
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FIG. 5: (Color dependence) Magnetic field dependence of x' 
and x" at T = 34 K in high magnetic fields. H* is shown for 
u)/2n = 1065 Hz. The inset shows fraction of grains in SSS, 
X, as a function of Ho/H^s- 



tion of the grains in SSS on the surface as a function of 
dc magnetic field, x{Ho) = J^'^ sm{e)Q{l.lHc2- Ha)de. 
The result of this calculation is shown in the inset to 
Fig. El The value H* /H^^ = 1.35 corresponds to the 
fraction of superconducting grains x « 0.66. One could 
infer that x" provides a more accurate onset for SSS, 
then H* — 12.5 kOe and x « 0.58. These values are in 
good agreement with the result for critical concentration 
of conducting clusters Xc = 0.6 in the 2d percolation site 
problem [2^ . We consider this as an additional argument 
in favor of the surface nature of the observed ac response. 

The adequateness of the GL model for the description 
of MgB2 was discussed recently in in frame of the Us- 
adel equations. In the dirty limit, the GL model is correct 
if the dimensionless parameter rj = Dc/2ttT is less than 



1. Here q « l/a^o is the average Fourier component of 
fs (x) in Eq. ||2Jl and Dc is the diffusion coefficient through 
which the tensor /t could be expressed 0. To be sure 
that calculations of Hc2 and Hc3 are correct we found 
parameter 77 for the most sensitive orientation of c-axis, 
Hq _L c and ip = tt/2 assuming that xq is the distance at 
which the order parameter in the 7r-band decreases by the 
factor of three. It was found that 7/ < 1 for T/Tc > 0.8 
and GL equations could be used at T/Tc > 0.8. An- 
other argument in favor of the GL equations is the very 
good agreement between the temperature dependence of 
Hcz/Hc2 for in ab plane calculated from Eq. (Q) and 
the result obtained in [22| on the basis of Usadel model. 

In summary, we have demonstrated the existence of 
SSS in polycrystalline MgB2. The SSS coexist with weak 
bulk magnetization. The ac response is linear at low am- 
plitudes of excitation field. The transition to SSS is of 
2d percolation character. The frequency dispersion for 
cij/27r in the interval of 5-1065 Hz shows a maximum in 
x", and x' is proportional to Inw. The superconduct- 
ing current in the grain depends on the orientation of 
the principal axes and is a function of the instant values 
of magnetic field and both k and kz- The relaxation of 
k and k^ to their equilibrium values with zero surface 
current determines the ac response of the grain. The 
subsequent averaging over all grains provides a complex 
character of the observed ac response. We found that GL 
equations adequately describe the equilibrium dc magne- 
tization of MgB2 samples at large field and could be used 
for analysis of SSS. It was shown that the ratio Hc3/Hc2 
weakly depends on the orientation of the c-axis whereas 
both Hc2 and Hc3 strongly depend on the orientation, 
which leads to the coexistence of SSS and weak bulk dc 
magnetic response. 
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